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DESCRIPTION OF THE INVENTION 

Cross-Reference to Related Application 

[001] This application claims the priority benefit of U.S. Provisional Application 
No. 60/259,121 , filed on December 29, 2000, which is incorporated herein by reference. 

Field of the Invention 

[002] The present invention relates to a method and apparatus for reducing 
interference in a received satellite signal using real-time signal processing. 

Background of the Invention 

[003] Navigational aides, such as those devices used in automobiles to assist 
drivers in locating destinations, have become very popular in recent years. These 
navigational aides work by receiving satellite signals from systems such as the Global 
Positioning System (GPS). GPS consists of 24 satellites that orbit the earth and 
transmit signals to these navigational aides. A navigational aide processes these 
signals to determine, for example, the location of a driver and, based on the driver's 
location, the navigational aides may provide directions to the driver's destination. In 
addition to navigational devices, GPS provides means for automatic vehicle location 
systems, aircraft landing systems, and precision timing systems. These devices have 
both commercial and military applications. 

[004] However, the satellite signals on which these devices rely are transmitted 
at a very low power level and are therefore susceptible to unintentional and intentional 
interference. Sources of unintentional interference include cellular phones and 
television stations transmitting antennas. Intentional interference (jamming) is 
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accomplished by intentionally producing signals to interfere with the satellite signals 
transmitted. 

[005] When interference occurs, the performance of devices that rely on the 
satellite signals degrades. To maintain or improve the performance of these devices in 
the presence of interference, GPS receivers must be designed to cancel or minimize the 
interference. 

[006] For a significant reduction of the effects of interference to a desired 
satellite signal, a hardware implementation of digital filters operating on analog-to-digital 
sampled data from the satellite receiver's intermediate frequency may be required. The 
digital filters require numerical coefficients that are derived from the incoming sampled 
data and are applied to the filters in real-time. However, in conventional systems the 
hardware needed to store the sampled data is very costly and computationally 
inefficient. Therefore, there is a need for a GPS receiver that may cancel or minimize 
interference in a cost effective and computationally simplified manner. 

SUMMARY OF THE INVENTION 

[007] A method for computing filter coefficients of a beamformer based on a 
segment of input samples. The method comprising the steps of dividing the segment of 
input samples into a plurality of blocks of input samples and receiving the plurality of 
blocks of input samples in a shared memory at a first rate wherein a first block of the 
plurality of blocks is received in a shared memory at a first time. The method further 
comprises reading the plurality of blocks of input samples from the shared memory at a 
second rate wherein the first block of the plurality of blocks is read from the shared 
memory at a second time. Still further, the method comprises computing a plurality of 
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partial covariance matrices for the plurality of blocks read from the shared memory and 
adding the plurality of partial covariance matrices. 

[008] Additional objects and advantages of the invention will be set forth in part 
in the description which follows, and in part will be obvious from the description, or may 
be learned by practice of the invention. The objects and advantages of the invention 
will be realized and attained by means of the elements and combinations particularly 
pointed out in the appended claims. 

[009] It is to be understood that both the foregoing general description and the 
following detailed description are exemplary and explanatory only and are not restrictive 
of the invention, as claimed. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[010] The accompanying drawings, which are incorporated in and constitute a 
part of this specification, illustrate several embodiments of the invention and together 
with the description, serve to explain the principles of the invention 

[01 1] Fig. 1 illustrates an exemplary embodiment of a spatial temporal adaptive 
processing (STAP) beamformer; 

[012] Fig. 2 illustrates an exemplary block diagram of a coefficient processor; 

[013] Fig. 3 illustrates an exemplary method for computing filter coefficients for 
the beamformer of Fig. 1 using the coefficient processor of Fig. 2; 

[014] Fig. 4 illustrates another exemplary block diagram of a coefficient 
processor; and 

[015] Fig. 5 illustrates an exemplary method for computing filter coefficients for 
the beamformer of Fig. 1 using the coefficient processor of Fig. 4. 
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DESCRIPTION OF THE EMBODIMENTS 

[016] Reference will now be made in detail to the present exemplary 
embodiments of the invention, examples of which are illustrated in the accompanying 
drawings. Wherever possible, the same reference numbers will be used throughout the 
drawings to refer to the same or like parts. 

[017] Fig. 1 depicts an exemplary embodiment of a spatial temporal adaptive 
processing (STAP) beamformer 100. The beamformer 100 may comprise N antenna 
elements 160 to receive one or more satellite signals 100A arriving at the antenna 
elements 160 in one or more directions. In addition to satellite signals 100A, the N 
antenna elements may also receive one or more interference signals 100B arriving at 
the antenna elements 160 from one or more directions. Each antenna element 160 is 
connected to a multiple tapped delay line structure 150 comprising M taps. Each of the 
multiple tapped delay line structures 150 may comprise M-1 delay elements 110, M 
multipliers 120, and an adder 130. Generally, each of the multiple tapped delay 
structures 150 may have an FIR structure. 

[018] Although not shown in Fig. 1, the signals 100A, 100B received by each 
antenna elements 160 may undergo preprocessing prior to being received by the 
multiple tapped delay line structures 150 and a coefficient processor 170. For example, 
signals received by each antenna elements 160 may be filtered by a preselection filter 
and down-converted to baseband or other IF frequency. Further, the signals may be 
filtered by a bandlimited filter and sampled by an analog-to-digital converter prior. Still 
further, the baseband filtered sampled signals may be further converted to complex 
baseband signals by digital demodulation or Hilbert transform type processing, for 
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example, prior to being input to the multiple tapped delay line structures 150 and the 
coefficient processor 170. 

[019] The beamformer 100 receives input signals 100A, 100B and computes 
filter coefficients, w nm , which are applied to the multiple tapped delay line structures 150 
for processing the input signals. The signals at the taps in the multiple tapped delay line 
structures 150 are weighted by corresponding filter coefficients, w nm , and the resulting 
products are summed by adders 130 producing an output signal for each multiple 
tapped delay line structure 150. The filter coefficients, w nm , are computed by the 
coefficient processor 170, which will be discussed in greater detail below. The outputs 
from each of the multiple tapped delay line structures 150 are then summed together by 
an adder 140 to generate output samples, y(k). 

[020] The output samples, y(k), of the beamformer 100 may be expressed by 
the following equation: 



N M 



[021] y(k) = YZ^n(k-{m-\)) 



(D 



n=\ m=l 



[022] where x n (k) denotes a complex input sample from the n-th antenna ! 
element at time k. It is assumed that at the n-th antenna element, the satellite signal is j 
multiplied by a factor e~ JK . The exponent factor, A n , depends on the angle of arrival of : 
the satellite signal, the carrier frequency of the satellite signal, and the position of the n- i 
th antenna element. More specifically, A n = o) c t n , where co c is the carrier frequency of j 
the satellite signal and r n is the inter-element time delay at antenna element n. If the i 
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steering vector for a given satellite direction is denoted by [e" M e" M ■ ■e~ jAN \, then the 
input samples, x„ (k) , may be expressed in the z-space by the following equation: 



[023] X n (z) = e- jA -e- jmT "V(z) 



(2) 



[024] where m is the baseband frequency of processed satellite signal received 
by the multiple tapped delay line structures 150 and a coefficient processor 170 and 
V(z) is the z-space representation of the satellite signal at a first antenna element. 

[025] To minimize the effects of interference signals 100B on the GPS receiver, 
the expected power, P, of the complex output samples, y(k), of the beamformer 100 
may be minimized according the following equation: 

V If M 
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[026] p = }= ^(fzi: - (» - i»Yzi:^(* - c/ -d) 



v «=1 m-l 
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[027] = E Vn(k -(m -lMk-(j -1)W 



(3) 
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[028] The expected power, p, may be simplified by rearranging the input 
samples and weights into the following (N x M) x 1 vectors: 
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[029] x(k) = 



x x {k) 
x l (k-l) 

Xl (4) 
x 2 (k) 

x 2 (4) 

x N (k) 
x N (k-l) 

%(4) 



(4) 
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[031] In matrix notation, the output samples, y{k), of the beamformer 100 may 
be expressed as follows: 

[032] y(k) = w H x(k) (6) 
[033] The resulting expected output power, P , is given by: 

[034] P =E{y(kf\=E{(w H x(k)lx H (k)w)} 
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[035] =w H E{x(k)x H (k)p =w H R~w 
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[036] where R~ is a covariance matrix. Before determining the minimum 

expected output power, P , at least one constraint may be imposed to avoid the trivial 
solution of zeros for the filter coefficients, w. Accordingly, equation (7) may be 
minimized subject to the following constraint: 

[037] C H w = F (8) 

[038] where C is a constraint weighting matrix and F is a constraint solution 
vector. The constraint weighting matrix, C, may be an (N x M) x L matrix and constraint 
solution vector, F , may be an L x 1 matrix. Accordingly, the interference minimization 
problem my be characterized as follows: 

[039] minimize P = w H R~w 

[040] subject to C H w = F 

[041] The filter coefficients, w, that may solve the interference minimization 
problem may be determined by the following equation: 

[042] w = R^c[c h R^cYf (9) 
[043] Conventional coefficient processors receive and store a predetermined ; 
number of input samples prior to computing the covariance matrix, R~ . For example, 
for a beamformer 100 having four antenna elements 160 (i.e., N = 4) and five taps (i.e., 
M = 5), the conventional coefficient processor may receive and store four thousand \ 
complex input samples prior to computing the covariance matrix, R~ . Because the 
input samples are complex, the coefficient processor may require a memory device 
having eight thousand memory locations. Further, conventional coefficient processors 
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are computationally intensive. In the above example, the conventional coefficient I 
processor computes the inversion of a 20 x 20 covariance matrix, R~ , and solves a 20- 

by-20 system of linear equations to compute the optimum filter coefficients, w . j 

! 

[044] Fig. 2 illustrates a first embodiment of the coefficient processor 170 ! 
according to the present invention. The coefficient processor 170 is a computationally ; 

i 

efficient processor for computing the covariance matrix, R~ . Further, the coefficient ! 

i 

processor 170 is less costly than conventional coefficient processors because j 
coefficient processor 170 requires less memory storage to compute the covariance \ 

i 

I 

matrix, R~ . Coefficient processor 170 may comprise a data input device 210 to receive ! 

j 

sampled input signals from the antenna elements 160, a shared memory device 220 to j 

store the input samples, a CPU 230 to compute the filter coefficients, w , a position input ; 

i 

device 240 to receive satellite position data, and an output device 250 to output the filter ■ 
coefficients, w, to the multiple tapped delay line structures 150 of the digital j 

l 

beamformer 100. The CPU 230 may compute the covariance matrix, R~ , by • 

i 

computing a plurality of partial covariance matrices, R~, and adding the partial I 

covariance matrices, RL , together to compute the covariance matrix, R~ . 

[045] More specifically, the elements, r^, of a covariance matrix, R~ , of an N x > 

S matrix where a sample in the N x S matrix is designated as ay, where the subscript i j 
refers to i-th row of the N x S matrix and j refers to the j-th column in the N x S matrix, ! 
may be determined based on the following equation: ! 



[046] r. =jy kl a kJ for l<i<S and 1 < j < S 



k=l 



(10) 
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[047] where the term a* ki is the complex conjugate of a . The computation of 

ki 

the elements, r.., of the covariance matrix, R~ , may be broken up into a sum of parts by 
dividing the N x S matrix into L sub-matrices where all L sub-matrices may have an M 
rows and S columns. Accordingly, the elements, r tj , of the covariance matrix, R~ , may 
be determined based on the following equation having L individual summation terms: 

M 2M LM 

[048] r ij =^a* ki a lg + £ a * kl a kJ +...+ ( 10 ) 



k=M+l 



k=(L-l)M+l 

[049] Each summation represents a partial covariance computation. 
Accordingly, the CPU 230 may compute the covariance matrix, R~ , by computing a 
plurality of partial covariance matrices, R'~ , and adding the partial covariance matrices, 
RL , together to compute the covariance matrix, R~ . The CPU 230 may then compute 
the filter coefficients, w , according to equation (9) above. 

[050] Fig. 3 illustrates a method used by the coefficient processor 170 of Fig. 2 
to compute the covariance matrix, R~ , for segments of input samples. Continuing with 
the example of a beamformer 100 having four antenna elements 160 (i.e., N = 4) and 
five taps (i.e., M = 5), the segments may consist of four thousand complex input 
samples grouped in eight blocks of 25 x 20 matrices. A block of 25 by 20 input samples 
may comprise twenty-five vectors having twenty input samples. The twenty input 
samples may consist of five input samples from each antenna element 160. More 
specifically, a block of 25 by 20 input samples may consist of the following input 
samples: 
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[051] 



jc, (* - 1 20) x, (* - 1 2 1) • • • x, (* - 1 24) x 2 (A - 120) • • • x 4 (k - 120) • • • x 4 (k - 124)' 

x,(£-5) x,(£-6) ■■■ Xl (k-9) x 2 (k-5) ••• x 4 (k-5) ■■■ x 4 (k-9) 
*,(*) x x {k-\) ... x,(*-4) x 2 (t) ••• x 4 (&) ••• x 4 (k-4) 



[052] As shown in Fig. 3, the shared memory device 220 and the CPU 230 
operate in parallel. That is, as the shared memory device 220 receives blocks of input 
samples for an i-th segment of input samples, the CPU 230 reads the blocks of input 
samples and processes the blocks of input samples to compute a covariance matrix, 
, for the i-th segment of input samples. As shown in Fig. 3, the shared memory 
device 220 may continuously receive input samples at a rate equal to 1125 
microseconds per block of input samples, for example. 

[053] The CPU 230 may not read a first block of input samples of the i-th 
segment to begin computing the covariance matrix, R~ , for the i-th segment until a time 
after the first block of input samples has been received by the shared memory device 
220. The purpose of the time delay may be to permit incoming input samples to fill the 
shared memory device 220 at a rate less than or equal to the rate the input samples are 
read by the CPU 230. In this way, blocks of input samples may always be available for 
processing by the CPU 230. To determine the number of blocks that may be stored in 
shared memory device 220 before the CPU 230 may read blocks of input data from the 
shared memory device 220 so that blocks of input samples may always be available for 
processing by the CPU 230, let the following variable be defined as follows: 

[054] N s = the number of blocks of samples in a segment of input samples; 

[055] N c = the number of blocks of samples in the shared memory device 220 
for a segment before a first block of input samples is read by the CPU 230; 
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[056] T w = the time to write a block of input samples to the shared memory 
device 220; 

[057] T r = the time to read a block of input samples from the shared memory 
device 220 by the CPU 230; and 

[058] T s = the time of one segment of input samples. 

[059] At a time t = N S T W , the number of blocks of input samples remaining in the 
shared memory device 230, N m , may be given by the following equation: 



[060] Nm = Ns - 



Tr 



(NsTw-NcTw) 



(11) 



[061] For blocks of input samples to be available for processing by the CPU : 
230, the number of blocks of input samples remaining in the shared memory device 
230, N m , at a time t = N S T W should be greater than or equal to zero. Accordingly, based 
on equation (1), the of number of blocks, N c , that may be stored in shared memory 
device 220 before the CPU 230 may read blocks of input data from the shared memory 
device 220 so that blocks of input samples may always be available for processing by , 
the CPU 230 may be determined based on the following equation: 



[062] Nc>M 



Tw 



(12) 



[063] For the embodiment of Fig. 3, where N s = 8, T w = 1 1 25 microseconds, T r = 
1100 microseconds, and T s = 10000 microseconds, N c may be at least 0.178 blocks. 
Accordingly, the CPU 230 may begin to read blocks of input data from the shared 
memory device 220 after 0.178 blocks of input samples have been received by the 
shared memory device 220 to ensure that blocks of input samples may always be 
available for processing by the CPU 230. If T w = 1250 microseconds and T r = 1100 
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microseconds, N c may be at least 0.960 blocks. If T w = T r , N c may be at least zero 
blocks. 

[064] Each time the CPU 230 reads a block of input samples from the shared 
memory device 220, it frees a memory cell in the shared memory device 220 to receive 
additional input samples from the data input device 210. In this way, the shared 
memory device 220 need not be capable of storing all four thousand complex samples 
needed to compute the covariance matrix, R~ . Accordingly, less memory storage is 
required to compute the covariance matrix, R~ . As illustrated in the examples above, 

the shared memory device 220 may only need to store a fraction of the total number of 
blocks of input signals for a segment. 

[065] Even if T w > T r , the shared memory device 220 may only need to store a 
fraction of the total number of blocks of input signals in a segment. For example, at i 
time t = N S T W and for N c = 0, the number of blocks of input samples remaining in the ! 
shared memory device 230, N m , may be given by the following equation: j 



[066] Nm = Ns 



1 — 



(13) 



[067] If T w = T r /2, the shared memory device 220 may only need to store one of 
the total blocks of input samples in a segment. 

[068] The CPU 230 computes the covariance matrix, R~ , for an i-th segment of 
input samples by computing a partial covariance matrix for each block of input samples 
for the i-th segment. The partial covariance matrix computation may be computed for 
the upper half of the covariance matrix, R~ , because the covariance matrix, R~ , is 
conjugate symmetric. Computing the partial covariance matrix in this way may result in 
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a saving of approximately half the computation time. When the CPU 230 has computed 
the partial covariance matrices for each block of input samples for the i-th segment of 
input samples, the CPU 230 adds the partial covariance matrices together to compute 
the upper half of the covariance matrix, R~ . The remaining half of the covariance 
matrix, R~ , is filled in by taking the complex conjugate of the i-th row and j-th column 
element and putting that value in the j-th row and i-th column of the covariance matrix, 
R~ . Once the covariance matrix, # H , is computed, the CPU device 230 may compute 
the filter coefficients, w, by first executing an LU decomposition algorithm which 
triangularizes the covariance matrix, R~ by decomposing it into the product of a lower- 
triangular matrix and an upper triangular matrix. Triangularization decomposes the 
linear system to be solved into two triangular systems of equations, which are solved 
recursively. Following LU decomposition, alternative algorithms may be used to 
compute the filter coefficients, w , for different operational conditions, w . For example, 

U.S. Application No. Serial No. , filed on even date herewith in the 

name of Khalil John Maalouf, Jeffrey Michael Ashe, and Naofal Al-Dhahir and entitled 
"A GPS Receiver Interference Nuller With No Satellite Signal Distortion," assigned to 
the assignee of the present application, which is hereby incorporated by reference, 
discloses algorithms that may be used to compute the filter coefficients, w . 

[069] Referring to Fig. 4, processing time for computing the covariance matrix, 
R^, may be further reduced by adding partial correlation processors 225 to the 
coefficient processor 170. Fig. 5 illustrates the process of the shared memory device 
220 operating in parallel with the partial correlation processors 225 and the CPU device 
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230. Continuing with the example of a beamformer 100 having four antenna elements 
160 (i.e., N = 4) and five taps (i.e., M = 5), the shared memory device 220 associates 
input samples in blocks of 25 by 20. The shared memory device 220 delivers at time t, 
an i-th block of 25 by 20 input samples to one of the partial correlation processors 225. 
The shared memory device 220 may continuously deliver blocks of 25 by 20 input 
samples to the partial correlation processors 225 evenly spaced in time every 375 
microsecond, for example. Each time the shared memory device 220 delivers a block 
of input samples to the CPU 230, it frees a memory cell in the shared memory device 
220 to receive additional input samples from the data input device 210. In this way, the 
shared memory device 220 need not be capable of storing all four thousand complex 
samples needed to compute the covariance matrix, R~ . 

[070] Each partial correlation processor 225 receives a block of 25 x 20 input 
samples from the shared memory device 220 and computes a partial covariance matrix. 
The number of partial correlation processors 225 may be chosen so that a partial 
correlation processor 225 is always available to begin processing a block of data 
received from the shared memory device 220. As shown in Fig. 5, four partial 
correlation processors 225 are provided. The first partial correlation processor 225 
receives the first block of input samples. While the first partial correlation processor 225 
processes the first block of input samples, the second partial correlation processor 225 
receives and processes the second block of input samples. While the first and second 
partial correlation processors 225 process the first and second block of input samples, 
the third processor 225 receives and processes the third block of input samples. 
Finally, while the first, second, and third processors 225 process the first, second, and 



15 



: ; 

5 q I 

1=5 



; - 

Li 



LAW OFFICES 

nnegan, Henderson, 
Farabow, Garrett, 

8 DUNNER, L.L.P. 

1300 I STREET, N. W. 
AS H I NGTONj DC 200 05 
202-4O8-4O00 



third block, the fourth processor 225 receives and processes the fourth block of input 
samples. When the fifth block of input samples is ready for processing, the first 
processor 225 is available to receive and process the fifth block. The remainder of the 
processors 225 are available to receive and process the remainder of the blocks of the 
segment of input samples. 

[071] When each partial correlation processor 225 finishes its computation, it 
delivers the result to the CPU device 230. When the CPU 230 receives the partial 
covariance matrices for a segment of input samples from the partial correlation 
processors 225, the CPU 230 adds the eight partial covariance matrices together to 
compute the covariance matrix, R~. The coefficient processor 170 of Fig. 5 may 
compute the covariance matrix, R~ , in a faster time than the coefficient processor 170 
of Fig. 3 because the multiple partial correlation processors 225 are available to read 
out and process the blocks of input samples quicker than the single CPU 230 of Fig. 3. 
For example, the coefficient processor 170 of Fig. 4 may reduce the processing time for 
computing the covariance matrix, R~ , to 3 milliseconds from 10 milliseconds using the 
coefficient processor of Fig. 2. 

[072] The coefficient processors 170 described above assumed a beamformer 
100 having four antenna elements 160 (i.e., N = 4) and five taps (i.e., M = 5). However, 
the coefficient processors 170 may be adapted for beamformers 100 having a variety of 
numbers of antenna elements 160 and taps. Adapting the coefficient processors 170 
for a variety of beamformers 100 will be obvious to those of ordinary skill in the art. 

[073] Other embodiments of the invention will be apparent to those skilled in 
the art from consideration of the specification and practice of the invention disclosed 
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herein. It is intended that the specification and examples be considered as exemplary 
only, with a true scope and spirit of the invention being indicated by the following claims. 
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